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. ABSTRACT 

The triple system HD 158926 (ASco) has been observed interferometrically with the 
\Q \ Sydney University Stellar Interferometer and the elements of the wide orbit have been 

determined. These are significantly more accurate than the previous elements found 
spectroscopically. The inclination of the wide orbit is consistent with the inclination 
previously found for the orbit of the close companion. The wide orbit also has low 
eccentricity, suggesting that the three stars were formed at the same time. 

The brightness ratio between the two B stars was also measured at A = 442 nm 
and 700 nm. The brightness ratio and colour index are consistent with the previous 
classification of A Sco A as B1.5 and A Sco B as B2. Evolutionary models show that the 
two stars lie on the main sequence. Since they have have the same age and luminosity 
class (IV) the mass-luminosity relation can be used to determine the mass ratio of the 
two stars: M B /M A = 0.76 ± 6.04. 

The spectroscopic data have been reanalyzed using the interferometric values for 
■ P, T, e and to, leading to revised values for aisini and the mass function. The 

individual masses can be found from the mass ratio, the mass function, spectrum 
synthesis and the requirement that the age of both components must be the same: 
M A = 10.4 ± 1.3 M Q and M B = 8.1 ± 1.0M@. 

The masses, angular semimajor axis and the period of the system can be used to 
determine the dynamical parallax. We find the distance to A Sco to be 112±5pc, which 
is approximately a factor of two closer than the HIPPARCOS value of 216 ± 42 pc. 

Key words: binaries: spectroscopic - binaries: visual - stars: fundamental parameters 
- techniques: interferometric 



1 INTRODUCTION 

The bright southern star A Scorpii (HD 158926, HR 6527, 
a 2 ooo = 17 h 33 m 36=52, <5 2 ooo = -37°06'13'.'8, V = 1.62) is 
classified as a B2IV+B single-lined spectroscopic binary sys- 
tem in the Bright Star Catalogue (Hofflcit & W arren 199111 . 
but it is in fact a triple system. Two recent papers by 
Uytterhoeven and colleagues l ] Uytterhoeven et al. 2 004a 
|Uytterhoeven et al. 2004b"l l present a detailed spectroscopic 

* E-mail: W.Tango@physics.usyd.edu.au 



study of the system based on fourteen years of data. We 
shall refer to these as Papers I & II. These papers also in- 
clude an extensive review of the literature on A Sco and here 
we present only a brief synopsis. 

The system comprises two B stars of similar mass that 
orbit each other with a period of approximately 1000 days 
(the 'wide' system) . The primary is a /3 Cephei type variable 
with a low mass companion that orbits the primary with 
a period of 6 days (the 'close' system). The spectroscopic 
elements for the two orbits are given in Paper I. The close 
binary system is eclipsing and in Paper I it is shown that this 
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constrains the inclination to the range 70° < i c io SC < 90°. 
An independent argument in Paper II based on frequency 
analysis also gives the same result. 

Paper II notes that the wide orbit cannot be accurately 
determined from the spectroscopy. We report here interfero- 
metric observations made with the Sydney University Stellar 
Interferometer (SUSI) on the wide system that results in a 
much improved orbit. 

In Section |21 we present the SUSI observations and the 
new interferometric orbit. 

In Section [3| the interferometric elements are used to 
determine a revised spectroscopic orbit. Interferometry also 
provides the brightness ratio between the distant component 
and the primary. The mass-luminosity relation can be used 
to find the mass ratio, and using the spectroscopically de- 
termined mass function, the individual masses can be calcu- 
lated. Having determined the masses the dynamical parallax 
for A Sco can be calculated. It differs significantly from the 
HIPPARCOS parallax. 

Following Papers I & II we use 'close' and 'wide' to 
distinguish the orbital elements associated with the 6-day 
and 1000-day periods seen in the spectrum of the primary. 
When no subscript is used the wide system is implied. 

In the case of the interferometric data we follow the 
scheme proposed in the Washington Multiplicity Catalog 
(WMC) Partkopf fc Mason 2003) ■ The primary will be de- 
noted by 'A', the close companion by 'a' and the wide com- 
ponent by 'B'. The notation a,Aa, for example, denotes the 
semimajor axis of the orbit of a relative to the primary while 
a.Aa,B is the semimajor axis of B relative to the centre of light 
of the close pair. 



2 THE SUSI OBSERVATIONS AND THE 
INTERFEROMETRIC ORBIT 

2.1 Observations with the 'blue' beam-combining 
optics 

Observations prior to 2004 were made at the blue wavelength 
of 442 nm using a bandwidth of 4nm. The optical layout 
was essentially the same as the one described in Davis et 
al. 1199911 . The data were reduced using the same procedure 
that was employed for the analysis of the binary system 
/3Cen lIDavis et al. 2005al . 

Baselines of 5, 10 and 20 metres were used. For these 
baselines the two B stars are essentially unresolved by SUSI. 
Because of atmospheric turbulence the fringe phase carries 
no useful information, and SUSI, like most optical interfer- 
ometers, measures a quantity proportional to the squared 
modulus of the complex degree of coherence 7: 

. |2 l+/3 2 +2/3cos(27tb-p/A) _ 

M iTTW (1) 

where b is the projected baseline, p is the vector separation 
between the two components of the binary as seen from the 
Earth and = Ig / 1 a ^ 1 is the brightness ratio. 

The SUSI data taken with the blue beam-combining 
system are 'uncalibrated'; i.e., no unresolved reference stars 
were observed to take into account the reduction in the ob- 
served fringe signal due to atmospheric and instrumental 
effects. The seeing conditions vary during the night, and ex- 
perience has shown that modelling this by a simple linear 



0.5 




0.0 1 — 1 — ' — 1 — ' — 1 — ' — 1 — ' — 1 — ' — 1 — ' — 1 — ' — 1 — ' — 1 — 1 
-5.0 -4.0 -3.0 -2.0 -1.0 0.0 1.0 2.0 3.0 4.0 

Hour Angle 



Figure 1. A typical set of data from one night (27 June 2000). 
Each point represents a measurement of \ V\ 2 . The data are fitted 
using equation J5J and the best fit is shown. The error bars denote 
the formal errors of the data points. The actual uncertainties are 
estimated to be 2.8 times larger. The point shown as an open 
circle was omitted from the final fit. 

time variation is usually adequate. We therefore assume that 
the observed signal will have the form 

|Vf = (C a ~C b t)\j\ 2 (2) 

where t is the hour angle of A Sco. The two empirical factors 
C a and Cb are free parameters in the model fitting proce- 
dure. The sign in equation (|5J reflects the fact that in general 
the seeing deteriorates during the night. 

It should be noted that the cosine term in equation 
depends on the position angle 9 and the separation p. 
When the raw \V\ 2 data are fitted using equation © the 
resultant estimates for p and 9 are not significantly affected 
by the fact that the data are uncalibrated. The estimate for 
the brightness ratio, on the other hand, may be subject to 
systematic errors. 

The presence of a cosine in equation (0 means that 
there is a 180° ambiguity in the position angle (a two- 
aperture interferometer cannot distinguish between p and 
— p). This has no impact on the interpretation of the results 
presented here. 

Equation Q does not take into account the presence 
of the close binary orbiting the primary. However, Paper I 
notes that the contribution of H a emission from weak T 
Tauri stars is at least 2 orders of magnitude fainter than the 
primary. Modulation due to the faint companion will be at 
most 1% and generally much smaller than this. This is less 
than the measurement noise and consequently the effect of 
the T Tauri component on the interferometric data will be 
negligible. 

The binary star measures were estimated from the \V\ 2 
data using essentially the same procedure described by Davis 
et al. (2005a). The data from each night were fitted using 
equation and the data for a typical night are shown 
m Fig. □ The formal error for each data point in general 
underestimates the actual scatter because the seeing varies 
throughout the night. Having found the best fit to the data 
the formal errors are rescaled to make \ 2 = v i where \ 2 
is the weighted sum of the squared residuals and v is the 
number of degrees of freedom. In effect, we assume that 
the actual uncertainties obey Gaussian statistics and all the 
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data points are similarly affected by seeing (in the exam- 
ple shown in Fig. 0the scaling factor was 2.8). Using the 
rescaled weighting factors Monte Carlo simulation is then 
used to determine the uncertainties in the fitting parame- 
ters. 

The assumption that all the points are equally affected 
by seeing is not strictly true, but if the seeing does not 
change dramatically during the night this assumption is, on 
average, realistic. The seeing is continually measured and 
logged, and points that are badly affected by seeing are ex- 
cluded from the final fit. 

A completely independent analysis was performed us- 
ing simulated annealing (Kirkpatr ick et al. 1983| l . Simulated 
annealing is more likely to converge on the true minimum of 
\ 2 rather than on a local minimum, unlike steepest descent 
methods. The simulated annealing results did not differ sig- 
nificantly from those obtained by our standard method and 
we are confident that the values of p, 8 and f3 and their un- 
certainties that are given in Table are the best fit to each 
night's data. The separation p and related quantities are all 
expressed in milliarcseconds (mas). 



2.2 Observations with the 'red' beam-combining 
optics 

The more recent SUSI data have been obtained with the 
red beam-combining optics. This system has been described 
in Davis et al. 12Q05.b| . The observations were made using 
baselines in the range 10-20 metres. Again the individual B 
stars are unresolved at these baselines and the effect of the 
six-day companion on the data will be negligible. 

The red observations were made at A = 700 nm using 
a bandwidth of 80 nm. With such a wide bandwidth the 
effects of "bandwidth smearing" cannot be ignored. Tango 
& Davis (2U02 j have discussed bandwidth smearing in the 
context of single stars and it is straightforward to generalize 
their results to arbitrary systems. In particular, if v = c/A 
is the optical frequency, S u is the spectral response of the 
interferometer and 



d x/„(x) 



(3) 



is the integrated monochromatic flux from the source (the 
integral is taken over the plane of the sky) , the broad band 
complex degree of coherence will be 



7 = 



duS v N u ~i[y) 
J °° duS v N v 



(4) 



where 7(f) is the quasimonochromatic complex degree of 
coherence defined by 
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NZ 



d xexp{— 27tib-xc V}/„(x) 



(•») 



The complex coherence 7(f) will vary with frequency (a) 
because of the frequency dependence of the Fourier kernel 
and (b) because I v (x) /N v is also a function of the frequency. 
In the case of a binary system with unresolved components 
this latter effect will be unimportant. 

We assume that the bandwidth is rectangular, having a 
width 5u centred on uq. Let Ao = c/v>o and make the usual 
approximation that 8u ~ c5A/Aq. For convenience define 



ip — 27tb-p/Ao- The equivalent of equation JQ including 
bandwidth smearing becomes: 

. |2 l + /3 2 sinc 2 (V><W2A„) + 2/3sinc(7MA/2A )cosi/> ,^ 

|T| = (TTW (6) 

In the case of A Sco the effect of bandwidth smearing can 
reduce |7| 2 by several percent. The effect is time-dependent, 
and is most significant when b and p are parallel. 

The red table observations were first processed using 
the "data pipeline" described in Davis et al. f2005bj. The 
output is a set of calibrated measures of |V| 2 . The visibility 
measures for each night were fitted using the same procedure 
adopted for the 'blue' observations. The only difference is 
that the red data are fitted using equation JSJ rather than 
equation The most noticeable difference between the 
blue and red data sets is that the errors in the red data 
points required very much less rescaling, confirming the fact 
that the use of calibrators significantly reduces the effects of 
seeing on the data. 

If the data are properly calibrated the maximum value 
of \V\ 2 , |V|maxj should be 1. On several nights, however, 
it was found that allowing VImax to be a free parameter 
gave significantly better fits, with | Vlmax < 1. The reduction 
cannot be explained by the partial resolution of the B stars 
or the modulation caused by the 6-day close companion. We 
believe that the effect is most likely instrumental but there 
are as yet not enough data to establish the exact cause. This 
effect has a negligible effect on the values of p and 8, but 
may affect f3. 

The results from the observations made with the red 
beam-combining system are given in Tabled 



2.3 The interferometric orbit 

The orbital elements were calculated, again using the same 
procedure that was employed for (5 Cen JDavis et al. 2005aTl . 
Each measure in Table^was given a weight according to the 
formula: W = (o 2 p +p 2 o 2 B )- 1/2 . The results are given in Table 
121 and the '\0 — C\' residuals are tabulated in Table Q 

The uncertainties in the orbital elements were estimated 
using Monte Carlo simulation using the best fit orbital el- 
ements and the weighting factors in Table to generate 
the simulated data sets. The resultant uncertainties are also 
listed in Tabled 

The period P, epoch of periastron T and eccentricity e 
differ significantly from those given in Paper I (see Tabic |3J, 
and we return to this in Section [J] 

The computed orbit and the observed measures are 
shown in Fig. [5] This figure also confirms the fact that any 
effects caused by the close companion must be smaller than 
the fitting errors. 

Our value of the inclination, iAa,B = 77?2±0°2, is com- 
pletely consistent with the range 70° ^ i c io SC ^ 90° given in 
Papers I & II and, together with the relatively low eccen- 
tricity of the interferometric orbit, supports the hypothesis 
that the B component did not join the system as a result of 
a tidal capture event. 



2.4 The brightness ratio and the mass ratio 

For the reasons discussed in Section l2.1l the brightness ra- 
tio f3 is subject to systematic errors associated with the 
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Table 1. The binary star measures determined with SUSI. The date, Julian date, the baseline used and the wavelength arc 
tabulated along with the measured values of p, 6, /3 and their uncertainties. The weighting factor W was used when fitting 
the orbit and is a measure of the accuracy of each point. The final column gives the residual between the observed data 
and the position calculated from the best fitting model. Using the WMC nomenclature, all these observations are classified 
as Aa,B. 



Date 


JD 


B 


A 


AA 


P 


<Tp 


( 


J 


°~6 




CT /3 


W 


\o-c\ 


CUT) 




(m) 


(nm) 


(nm) 


(mas) 


(mas) 
















(mas ) 


(mas) 


1 QQQ 0^ 1 ^ 


94^1 ^14 1 


10 


442 


4 


45.86 


0.19 


265 


?92 





?06 


0.532 





.008 


5 


0.51 


1 QQQ 0^17 


94^1 ^16 1 

iiTiUlulU. -L 


5 


442 


4 


45.66 


0.27 


266 


?37 





?08 


0.596 


0. 


.009 


3 6 


0.44 


1 QQQ 06 04 


94^1 ^34 1 


5 


442 


4 


46.12 


0.18 


267 


?95 





?07 


0.649 


0. 


.009 


5 2 


0.73 


1 QQQ 06 99 


94^1 3^9 1 


5 


442 


4 


46.79 


0.97 


268 


?93 





?27 


0.540 





.013 


1 


0.50 


1 QQQ 06 9^ 
iyyy uu 6o 


94^1 3^ 1 


5 


442 


4 


46.25 


0.83 


268 


?50 





?26 


0.580 


0. 


.023 


1 2 


1.11 


1 QQQ OK 1 9 
J. yyy-uo ± ^ 


94^1 40^ 1 


5 


442 


4 


44.53 


1.82 


273 


?09 





"13 


0.482 


0. 


.020 


5 


1.40 


£UUU"UJ" J. -L 


94^1 676 1 

£<-tO J.U 1 U.l 


10 


442 


4 


25.73 


0.21 


72 


?46 





?19 


0.553 


0. 


.015 


4 5 


0.34 


9000 0^19 


94^1 677 1 

Z^iU J_U 1 1 .1 


10 


442 


4 


26.08 


0.14 


72 


?73 





?10 


0.596 





.006 


6 6 


0.29 


9000 0^1^ 
iUUU uo ±0 


94^1 680 1 


10 


442 


4 


26.96 


0.16 


73 


?95 





"11 


0.590 





.019 


5 9 


0.16 


9000 O^i 1 R 
^UUU UJ ±o 


94^1 68^ 1 


10 


442 


4 


28.33 


0.11 


74 


?74 





?09 


0.609 


0. 


.003 


8.4 


0.52 


9000 O^i 97 

£UUU"Uul"« ( 


94^11 6Q9 1 


10 


442 


4 


30.23 


0.23 


76 


?49 





?19 


0.537 


0. 


010 


4 o 


0.18 


9000-06-09 


2451698 1 


10 


442 


4 


31.43 


0.14 


77 


?56 





?13 


0.528 





.005 


6.4 


0.06 


9000 06 0^ 
iUUU UU uo 


94^1 6QQ 1 


10 


442 


4 


31.45 


0.10 


77 


?57 





?07 


0.600 


0. 


.002 


9 


0.29 


9000-06-04 

iUUU UU Urt 


94^1 700 1 

ti'-i.xj J_ 1 UU. J- 


10 


442 


4 


31.25 


0.10 


77 


?72 





" 05 


0.575 


0. 


.005 


10.0 


0.72 


2000-06-07 


2451703.1 


10 


442 


4 


32.41 


0.13 


78 


?35 





?06 


0.529 


0. 


.008 


7.3 


0.26 


2000-06-11 


2451707.1 


10 


442 


4 


33.45 


0.28 


78 


?83 





?T4 


0.468 


0. 


.008 


3.4 


0.21 


2000-06-27 


2451723.1 


10 


442 


4 


36.94 


0.14 


81 


?39 





?07 


0.613 





.014 


6.7 


0.10 


2000-06-30 


2451726.1 


10 


442 


4 


37.44 


0.24 


81 


?83 





?14 


0.455 


0. 


010 


3.9 


0.20 


2000-07-01 


2451727.1 


5 


442 


4 


40.00 


0.40 


82 


?08 





?18 


0.588 


0. 


.014 


2.4 


2.16 


2000-07-22 


2451748.1 


5 


442 


4 


43.54 


0.21 


84 


?10 





?09 


0.560 





.009 


4.4 


1.83 


2000-08-02 


2451759.1 


5 


442 


4 


44.30 


0.26 


85 


?30 





?16 


0.543 


0. 


.012 


3.5 


0.83 


2001-06-22 


2452083.1 


10 


442 


4 


18.80 


0.17 


129 


?16 





"13 


0.534 


0. 


.012 


4.5 


1.14 


2001-08-10 


2452132.1 


10 


442 


4 


12.97 


0.16 


163 


?25 


1 


?61 


0.563 


0. 


.018 


2.5 


0.62 


2001-08-19 


2452141.1 


20 


442 


4 


12.58 


0.04 


174 


?46 





?43 


0.476 





.009 


9.6 


0.54 


2004-04-15 


2453111.2 


20 


700 


80 


24.05 


0.20 


118 


?92 





?24 


0.544 





.001 


4.5 


0.69 


2004-04-20 


2453116.2 


20 


700 


80 


23.38 


0.24 


119 


?37 





?51 


0.518 a 





.028 


3.1 


0.34 


2004-06-24 


2453181.0 


20 


700 


80 


12.77 


0.03 


160 


?40 





?30 


0.488 


0. 


.006 


13.3 


0.26 


2005-04-17 


2453478.2 


20 


700 


80 


47.42 


0.20 


270 


?36 





?05 


0.516 


0. 


.014 


4.8 


0.18 


2005-07-18 


2453570.0 


15 


700 


80 


37.83 


0.37 


279 


?45 





?10 


0.635 a 





.019 


2.6 


0.17 


2005-07-21 


2453573.0 


15 


700 


80 


37.18 


0.23 


279 


?76 





?05 


0.536 a 





.013 


4.2 


0.24 


2005-08-17 


2453600.0 


15 


700 


80 


32.23 


0.28 


283 


?68 





?07 


0.543° 





.008 


3.5 


0.45 



a Knax < 1 f° r these sets of data. The other observations at 700 nm were fitted using \ V\ 2 = 1. This is discussed in Section 
l2~2l 



of approximately 0.023 magnitudes and a period of 5^1. 
UShobbrook fc Lomb 19721 IShobbrook fc Lomb 19751 . 

Shobbrook & Lomb 1197511 also observed periodic vari- 
ations in colour. The amplitudes of the variation were 
(&-«) = 0.002 ± 0.001 and [u-b) = 0.0098 ± 0.0012. 

The measured brightness ratios and their uncertainties 
are plotted in Fig. |3 The data fall into two groups: the 
earlier set of points were observed using A = 442 nm while 
the more recent set was obtained using A = 700 nm. The 
two lines represent the approximate spread in j3 that would 
be expected from the intrinsic variability (~ ±0.012 mag). 
The error bars indicate the statistical uncertainties for each 
observation, but in the present context they are misleading, 
since there are significant night-to-night variations due to 
the intrinsic variability of A Sco as well as the instrumental 
effects (calibration errors) previously discussed. In our view 
the average and standard error of the blue and red data sets 
represent the best estimates of the brightness ratios. 

The means and standard deviations for the brightness 
ratio for each set, and the corresponding magnitude differ- 



Table 2. The orbital elements of A Sco determined from the SUSI 
observations. All the elements refer to the Aa,B orbit. 



Element Value, uncertainty and units 
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1052.8 ± 1.2 
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51562.3 ± 2.8 
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See note. 
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(a"f/P 2 


(108 ± 2) x 10" 12 


asec 3 d -2 



Note: This is the argument of periastron for the orbit of the B 
component relative to the primary. The corresponding quantity 
in Table |3 refers to the orbit of the primary around the system 
barycentre. The two angles differ by 180°. 

way the \V\ 2 data are normalized. As well, the system is 
a f3 Cephei variable with a photometric amplitude in v 
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Figure 2. The orbit of A Sco. The scale is in millarcsec and the orientation conforms to standard practice (Hcintz 19781. Different 
symbols are used to denote successive orbits, referred to the epoch of periastron T. The circles indicate data taken in 1999. The triangles 
are the data from 2000 and 2001. The diamonds denote points on the next orbit during 2004 and 2005. Open symbols are used for 
observations at 442 nm; filled symbols denote 700 nm observations. 
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Figure 3. The brightness ratios, taken from Tabled The two 
lines denote the range in that is expected due to the intrinsic 
variability of the system. The data prior to JD 2 452 500 were 
measured at A = 442 nm. The more recent data were measured 
at 700 nm. 

ences, are 

/3 44 2 =0.55 ±0.05 Am 44 2 = 0.65 ±0.10 (7) 

/3 700 = 0.54 ± 0.05 Am 7 oo = 0.67 ± 0.10 (8) 

The difference in the colour index (771442 — m,7oo) between A 
Sco A and B is —0.02 ±0.14. Both the magnitude and colour 
index differences are consistent with the relative spectral 
classification of A Sco A as B1.5 IV and A Sco B as B2 V 
given in Paper I. 

Although the A and B components have apparently dif- 
ferent luminosities, due to the complexity of this system and 
the fast rotation there is an uncertainty of up to 0.2 dex in 
log<?, making it difficult to establish the luminosity class 
with certainty. However, the inclination of the long period 
orbit suggests that the three stars formed together and A Sco 
A and B are thus both on the main sequence. The bright- 
ness ratio then implies that the effective temperatures can 



be fixed at 25 000 K and 21 000 K. The contribution of the 
T Tauri star to the luminosity of the primary is negligible 
and the empirical mass- luminosity relation log(Ls / La) = 
(3.51±0.14) log(M s /M A ) HGriffiths et al. 198811 can be used 
to find the mass ratio for A Sco A and B. The ratio of 
the luminosities has been obtained from the V magnitude 
difference, taken to be 0.66±0.10 by interpolation between 
Am 4 42 and A771700, and the difference in bolometric correc- 
tion (BC). The latter has been obtained from interpolation 
in the tables by Flower (I1996t for main-sequence B stars. 
Assuming an uncertainty of ±1000 K in the effective tem- 
peratures as obtained in Paper I, we arrive at (BCa - BCs) 
= -0.40±0.10. Thus the mass ratio for A Sco A and B is: 

Mr 

^-=0.76 ±0.04 (9) 

M A v ' 

This is consistent with the value of Mb /Ma = 0.84 ± 0.12 
estimated from the spectrum synthesis analysis presented in 
Paper I. 



3 THE COMBINED INTERFEROMETRIC 
AND SPECTROSCOPIC ORBIT 

The spectroscopic orbital elements were determined in Pa- 
per I using two different codes, FOTEL and VCURVE. The 
upper panel of Fig.[I]prots the RV data against the FOTEL 
phase. The solid line is the calculated RV using the FOTEL 
elements (this panel is the same as the lower panel in Fig. 2 
of Paper I). 

The interferometric values for P, T, e and ui are much 
more accurate than the spectroscopically determined ones, 
and the original RV data were refitted using these values, 
giving ifwidc = 22.8 ± 0.4km-s -1 . Table [3] tabulates the ele- 
ments given in Paper I and our revised values, and the RV 
amplitude calculated using these elements is shown in the 
lower panel of Fig. 0] 
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Figure 4. The radial velocity data for the 'wide' system in A Sco. 
The RV due to the 'close' system has been removed. In the upper 
panel the data are plotted against the FOTEL phase and the RV 
curve is calculated from the FOTEL elements in Table |3J. The 
same RV data are plotted in the lower panel against the SUSI 
phase, with the RV curve calculated using the elements given in 
the final column of Table . 



Table 3. The second column lists the elements found in Paper I 
using the FOTEL code. The third column tabulates the elements 
used in the final fit: P, T, e and lu are based on the SUSI; K 
has been recomputed using these values. The uncertainties are 
the formal uncertainties of the fit and do not include systematic 
effects. 



Element 


FOTEL 


FOTEL+SUSI 


P(d) 


1082 ± 3 


1052.8 ±1.2 


K (km-s" 1 ) 


24.7 ±0.4 


22.8 ±0.4 


T (MJD) 


51 731.5 ±29 


51 562.3 ±2.8 


e 


0.23 ±0.03 


0.121 ±0.005 




311° ±11° 


254?8±0?9 



Some words of caution on the accuracy of the RV data 
are in order. The use of two types of spectra of very dif- 
ferent nature and quality forced the authors of Paper I to 
use different methodology to derive the RV values. In doing 
so, they ignored the fact that the B component contributes 
to the spectral lines used for the RV computation, given 
that the binary is observed to be single-lined. The light, and 
hence mass ratio, found from the interferometry is an or- 
der of magnitude more precise than the one in Paper I and 
thus allows us to place better constraints on the individual 
masses of the A and B components. The ranges listed in Pa- 
per I must be revised, as they do not include the systematic 



Table 4. System parameters that depend on the RV amplitude. 
The lower limit for if w id c is taken from Table 151 The upper limit 
is based on our best estimates of the effects of the wide component 
on the RV data and from the interferometric data. 



Element 


Range 






[22.4, 25.7] 


km-s 


ai sini widc 


[2.15, 2.47] 


AU 


ai 


[2.21, 2.53] 


AU 


f(M) 


[1.20, 1.82] 


M Q 


M A 


[9.39, 14.25] 


M 


M b 


[7.05, 10.69] 


M 



effects due to the B component on the spectral lines used to 
determine the amplitude of the RV. 

To make the best use of the interferometric results it is 
important to understand the uncertainties and possible sys- 
tematic effects in the RV data. In general, the contribution 
of the B component to the lines leads to an underestimation 
of the RV of the primary. As a consequence, the derived 
quantities asini and the mass function are also underesti- 
mated. This fact was clearly demonstrated and corrected 
for by spectral disentangling in the case of the double-lined 
binary /3Cen (Ausseloos 2006} , where it implied a 10% in- 
crease in the individual component masses. We are unable 
to apply such a treatment to A Sco because the disentangling 
fails in this single-lined case with three stellar components, 
of which at least one is oscillating (see Paper II) . Neverthe- 
less, we can estimate the systematic errors of the RVs for 
the two sets of data from which they were derived: 

(i) the high S/N single-order spectra including the Si ill 
4553A line taken with CAT/CES; 

(ii) the low S/N echelle spectra taken with Eu- 
ler/CORALIE. 

The resolution for both sets of data is nearly the same (see 
Paper I for details). 

From set (i) we derive that the lines due to the B com- 
ponent are completely situated within those of the primary 
(see Paper I and Waelkens 0,990 ) for appropriate plots of 
the lines). Moreover, from the brightness ratio and the tem- 
perature ranges listed in Paper I we derive that the A Sco B 
contributes 40% to the total equivalent width of the Si ill line 
used to determine the RVs. This places strong constraints 
on the upper limit to K wl d c : the fact that the B component 
is not seen in the wings of the high-quality CAT profiles 
implies that K^-^e < 25.7 km-s -1 (which we deduce from 
merging the contributions of the two stars to the Sim line 
with the appropriate EW ratios and v sin i estimates of Pa- 
pers I and II) . This suggests a conservative range for i\wid c is 
[22.4, 25.7] km-s -1 . Table[I]tabulates the system parameters 
that are affected by the uncertainties in the radial veloc- 
ity amplitude. The accurate value of the mass ratio given 
in equation © has been used to calculate the individual 
masses from the mass function. 

The mass ranges listed in Tabled] although compatible 
with those given in Paper I, are very broad. We tried to con- 
strain them further by taking into account the results from 
the spectrum synthesis presented in Paper I, in which the 
contributions from the individual components were properly 
taken into account to fit H, He, and Si lines in data set (ii). 
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Table 5. Summary of the stellar parameters of the three com- 
pocnts of A Sco. The mass of the a component is taken from Pa- 
per I; the masses of A Sco A and B are based on the work presented 
here (compare with Table 4 of Paper I). 





Mass (M©) 


T cft (K) 


log 3 




A 


10.4 ± 1.3 


25 000 ± 1000 


3.8 ±0.1 


B1.5IV 


a 


1.8 ±0.2 






pre-MS 


B 


8.1 ± 1.0 


21 000 ± 1 000 


4.0 ±0.1 


B2IV 



We did this by scanning the data base of standard stellar 
evolution models computed by Ausseloos et al. ( 2003J. We 
selected models from the ZAMS to the TAMS with X = 0.70 
and Z in the range 0.012 to 0.030. Models with the three 
values 0.0, 0.1 and 0.2 for the core overshooting parameter 
(expressed in units of the local pressure scale height) were 
included. While scanning the database, we required that the 
two stars not only satisfy the mass ranges given in Table |1J 
but also that they have the appropriate luminosity ratio 
given by the interferometry, and have logT e g and loggr in 
the ranges found from the spectrum synthesis done in Pa- 
per I. Moreover, we required that their ages be the same to 
within 0.1% (following the interferometric result that the or- 
bital inclinations of the A and B components are consistent 
with each other). In this way, we find the acceptable range 
for the mass of the primary Ma to be [9.0, 11.7] Mq and the 
range for Mb to be [7.0,9.1] M . 

All the acceptable model combinations at the lower 
mass ends of these intervals have an age below 10 7 yr cor- 
responding to less than 60% of the main sequence lifetime, 
while those at the upper mass end correspond to an age be- 
low 13 x 10 r yr which is less than 30% of the main sequence 
duration. For the lower mass end, this is consistent with the 
suggestion made in Paper I that A Sco contains a pre-MS 
star. Indeed, from Palla & Stahler (19933 we n °te that a pre- 
MS star with a mass of ~ 1.5 Mq remains about 10 7 yr in 
its pre-MS stage when a B-type star formed from the same 
accretion disk has reached the ZAMS. Table summarizes 
our conclusions regarding the three components of A Sco. 



3.1 The dynamical parallax 

The dynamical parallax is 



365.25 



2/3 



(M Aa + M B y/3 



(a") 3 

p2 



1/3 



(10) 



where Ma„ is the combined mass of the primary (A) and 
the close companion (a) and Mb the mass of the B com- 
ponent. We have used the estimates given in Table |5] and 
the interferometric value of (a") 3 /P 2 (Table [5| to find the 
distance to A Sco: 



D = l/ir d = 112 ±5 pc 



(11) 



Our distance is smaller than the HIPPARCOS value 
of 216 ± 42 pc UPerryman et al. 1997^ by nearly a factor of 
two (i.e., 2.5 times the HIPPARCOS standard deviation). 
It has been noted (dc Zccu w^t al. 19991 that the motion 
in binary systems can significantly bias HIPPARCOS par- 
allaxes and a similar phenomenon has been observed with 
/3Cen lIDavis et al. 2005al . 



The revised parallax may have implications regarding 
the membership of A Sco in the Scorpius-Centaurus-Lupus- 
Crux complex of OB associations. Brown & Verscheuren 
( 1997 ) suggested that it was a member of the Upper Scorpius 
(US) association, but de Zeeuw et al. ©991 did not include 
A Sco as a secure member of the US association based on 
an analysis of HIPPARCOS positions, parallaxes and proper 
motions. The revised parallax corresponds to the Lower Cen- 
taurus Crux (LCC) association, which lies at a distance of 
116 ± 2 pc. The galactic coordinates of A Sco, however, still 
exclude it as a definite member of the association. 



4 SUMMARY 

Observations with SUSI have been used to establish an inter- 
ferometric orbit for A Sco. The accuracy of the period, epoch 
of periastron, eccentricity and argument of periastron are 
significantly better than those previously found (Paper I). 
We have combined the interferometric and spectroscopic el- 
ements to provide an orbital solution which is consistent 
with both the interferometric and spectroscopic data. The 
inclination of the orbit for the wide component is consistent 
with that previously found for the close companion, and the 
eccentricity of its orbit is quite low. This is strong evidence 
that the system was formed from a common accretion disk 
rather than through tidal capture. 

The brightness ratio allows the mass ratio of the two B 
stars to be estimated from the mass-luminosity relation. The 
estimated mass ratio of 0.76 ± 0.04 agrees with the masses 
determined by the synthetic spectrum analysis presented in 
Paper I. The limiting factor in the determination of the in- 
dividual masses and the semimajor axis is clearly the large 
uncertainties of the RV values. This is due to our inability 
to disentangle the contributions of the primary and wide 
companion to the lines in the overall spectrum. While so- 
phisticated disentangling methodology is available in the lit- 
erature, e.g. Hadrava (19951, its application to single-lined 
binaries has not yet been achieved to our knowledge. In the 
case of A Sco, the j3 Cep oscillations are an additional factor 
that complicates the disentangling. Nevertheless, we have 
taken a conservative approach to estimate the mass ranges 
from the spectroscopic orbital RV and we have constrained 
them further from spectrum synthesis results and from re- 
quiring equal ages and the appropriate brightness ratio from 
the SUSI data. In this way, we obtained a relative precision 
of about 12% for the masses of both A Sco A and B. Masses 
of B-type stars are not generally known to such a precision. 

Accurate mass estimates of pulsating massive stars are 
important for asteroseismic investigations of these stars. In- 
deed, the oscillation frequencies of well identified modes in 
principle allow the fine-tuning of the physics of the evolu- 
tionary models (see Aerts et al. 1200 3 1 for an example). For 
this to be effective when only a few oscillation modes have 
been identified we need an independent high-precision esti- 
mate of the effective temperature and of the mass of the star. 
This is generally not available for single B stars. The A Sco 
Aa system is one of the bright close binaries with a /3Cep 
component and it is only the second one, after /3Cen, for 
which a full coverage of the orbit with interferometric data is 
available. The results presented in this paper therefore con- 
stitute a fruitful starting point for seismic modelling of the 
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primary of A Sco, particularly now that the WIRE satellite 
has provided several more oscillation frequencies compared 
to those found in Paper II ijBruntt fc Buzasi 2005|l 

The mass estimates, when combined with the period 
and angular semimajor axis found from the interferometric 
orbit, allow the dynamical parallax to be calculated. We find 
the distance to the A Sco system of 112 ± 5 pc, which is al- 
most a factor of two smaller than the HIPPARCOS distance 
of 216±42 pc. The reason for the large discrepancy is almost 
certainly due to the B component since it is of comparable 
brightness to the primary and the orbital period is approxi- 
mately three years. The fact that many of the B stars in the 
Sco-Cen association are multiple systems suggests that the 
HIPPARCOS parallaxes for these stars must be used with 
care. 
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